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Abstract 

This short note presents a peculiar generalization of the Riemann hy- 
pothesis, as the action of the permutation group on the elements of con- 
tinued fractions. The problem is difficult to attack through traditional 
analytic techniques, and thus this note focuses on providing a numerical 
survey. These results indicate a broad class of previously unexamined 
functions may obey the Riemann hypothesis in general, and even share 
the non-trivial zeros in particular. 



1 Introduction 

The Riemann zeta can be expressed as an integral in the following formfT]: 
C(s) - — s f hix^x'-^dx 

S - 1 ./n 



where 

h[x) = -- 



X 



is sometimes called the Gauss map. Here, [yj is the largest integer less than or 
equal to y. The Gauss map lops of a digit in the continued fraction expansion 
of X. If one writes the continued fraction expansion^] of < a; < 1 as 



[01,02,03, 



Ol 



"3 + - 



then h{x) = [02, 03, • • • ]. It is a (reverse) shift operator on the continued fraction 
expansion. The study of continued fractions is interesting because of their prox- 
imity to various problems of fractals and dynamical systems, including ties to 
Farey fractions, the modular group via the Cantor set, and the Minkowski ques- 
tion mark function. More immediately, the transfer (Ruelle-Frobenius-Perron) 
operator of the Gauss map is the Gauss-Kuzmin-Wirsing operator, which has a 
number of interesting properties of its own. 

In this note, we observe that the Gauss map can also be thought of as one 
particular element of the permutation group acting on an infinite dimensional 
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representation of the real numbers. Thus, we are lead to contemplate other 
elements of the permutation group. Specifically, consider the permutation op- 
erator 

Sp,q : [fll, ■ • ■ , Op, ■ • ■ , flq, • • • ] i-> [oi, ■ ■ ■ ,aq, ■ ■ ■ ,ap, ■ ■ ■] (1) 

which exchanges the p'th and g'th digits of the continued fraction expansion of x. 
Note that Sp,q maps the unit interval onto the unit interval, and is discontinuous 
at a countably infinite number of points. The generalized Riemann hypothesis 
that we choose to explore concerns the zeros of the integral 

Cp.,(s) = - Sp,g{x)x''^dx (2) 

s ~ i Jo 

That is, is it possible that the (non-trivial) zeros of Cp,q{s) all lie on the dis — ^ 
axis? This text is devoted to the exploration of this possibility. 

This integral can be evaluated, but only with considerable difficulty. The 
function Sp^q (x) is a piece-wise combination of unimodular Mobius functions, 
that is, a piece- wise collection of 

ax + b 
cx + d 

with integer values of a, &, c and d, such that ad— be — ±1. This structure follows 
from the general theory of continued fractions Section |3] below provides a 
set of graphs visualizing Sp^q{x), showing their piece- wise nature. The functions 
are clearly "fractal" or self-similar; a general discussion of the self-similarity is 
given in [1]. Because of the piece- wise structure, the integral can be evaluated 
analytically. The result is a messy set of nested sums that don't seem to provide 
any particular insight. One might expect that these sums would simplify, or at 
least, speed up numeric evaluation, but they don't even seem very useful for 
that. Details are given in section O 

It is easiest, at first, to perform the integration numerically. Treating the in- 
tegral as a very simple Newtonian integration sum with equally spaced abscissas, 
one gets: 

CN;p,q (■^) = 7^ - ^ E ^P^. 

This summation converges very slowly and very noisily as oo but a nu- 

merical evaluation for large TV can give some basic insight. First and foremost, 
one discovers that, at least qualitatively, the sums CN-p.q (s) behave very much 
like the corresponding sum for the Riemann zeta. This provides the simplest 
evidence in support of the hypothesis. A graphical visualization is provided in 
section |3l 

2 Continued Fractions; Permutation Groups 

Continued fractions provide a representation of the real numbers in the infinite 
Cartesian product space = M x Mx... where M = N U {oo} and N = Z+ is 
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the set of positive integers. The continued fraction brackets are a function that 
maps this product space to the real numbers, i.e. [] : N"— ?> R. The mapping 
is, strictly-speaking, surjective, yet, in a certain sense, is "almost everywhere" 
bijective. That is, every non-rational real number corresponds to a unique 
continued fraction; only the rationals have multiple expansions. 

The above definition of a continued fraction differs slightly from the conven- 
tional norm. Usually, continued fractions are considered to be of finite length 
for the rationals, and infinite-length for the non-rational reals. This can create 
confusion when discussing Sp^q{x) when either p or g are larger than the length 
of the continued fraction. The trick of adjoining infinity to the naturals is the 
most straight-forward way of avoiding this problem. The cost of this trick is 
that a "very large part" of the Cartesian product maps to the rationals. Thus, 
for example, the continued fraction expansion for is = [oo, 02, 03, ...] for any 
arbitrary positive integers a^. All rationals suffer in the same way: a finite 
continued fraction is just one where 00 appears somewhere in the expansion, 
i.e. [ai, 02, Oat] = [oi, 02, oat, cxa, aAr+2, ...]. This can be papered over by 
considering the reals as the quotient space of sequences modulo the kernel of [ ] . 
That is, define a quotient space § = M'^/ker[], so that the operator [] : § — > M 
is 1-1 and onto in this quotient space. Since this map is bijective, it is invert- 
ible, and thus, any function on the reals R and be expressed as an equivalent 
function on S, and vice-versa: functions on the sequence space M'^ correspond 
to functions on the reals. 

Elements of the (infinite-dimensional) permutation group acting on S are 
given by equation [1] The Sp^q{x) are discontinuous, but only on a countable 
set of points. That is, for fixed p, g, the function Sp^q{x) is discontinuous on 
the proverbial "set of measure zero", and so one can easily form well-defined 
integrals using ordinary techniques. 

3 Structural overview 

A graphical presentation of the functions 5i.2 (a;), 5i.3 (a;), 6*1,4 {x) and 52.3 {x) 
are shown in figure[TJ The fractal, self-similar nature of these functions is readily 
apparent. 

Before exploring the sums[3]in detail, it is worth examining the corresponding 
sum for the Riemann zeta, so as to have a baseline to compare to. Define the 
Riemann sum as 

As a function of N , this sum converges very slowly to the Riemann zeta 
function. As N increases, the traditional zeros become visible, although, or 
finite N , they do not sit exactly on the 5Rs = | critical line. The traditional 
zeros are also accompanied by a larger number of "artifact zeros": closely-spaced 
zeros, with a much smaller residue, whose spacing (and residue) decreases as N 
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Figure 1: Continued Fraction Permutations 
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The above figures show a set of graphs for various continued fraction per- 
mutations Sp^q{x) for the swap operator Sp^q : [ai,--- ,ap,--- ,aq,---] i— J> 
[ai, • • • , flq, ■ • • , flp, • • • ]. A countable number of discontinuities are easily seen. 
It should also be clear that the density of these discontinuities make it quite dif- 
ficult to work with these functions using simple-minded numerical approaches. 
These figures posses a fractal self-similarity; this self-similarity is presumably 
key to obtaining workable analytic expressions. 
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rises. Visually, these are very easy to distinguish; the figure [2] show both the 
traditional and the artifact zeros clearly. 

The convergence properties of this sum, as well as the nature of the "artifact 
zeros", can be understood by exploring the phase of along the critical axis. 
This is shown in figure [31 Truncation at finite N introduces a quasi-sinusoidal 
perturbation overlaid on the "true" Riemann zeta function. 

Armed with this overview, one may now explore the sum Cn;1,2- As figure^] 
shows, the simplest numerical evidence seems to directly support the hypothesis. 
The evidence for the hypothesis, in the case of CAf;i,3 s-nd CAf;2,3 is less clear, as 
seen in figure [S] On the other hand, the numeric sums, for any fixed N, are 
far less "accurate" than that for Cn;1.2- That is, the discontinuities for 81^3 and 
'S'2,3 are very finely spaced, and the summation does a considerably poorer job 
of capturing these. 



4 A Remarkable Shadow 

The numerical evaluation of sums requires a certain amount of cross-checking 
so as to catch errors. A particularly curious sum to evaluate is 



ri=0 ^ 



which does no permutation at all. That is, one effectively has ?7Ar(s) = CN:p.p{s) 
for any integer p. The limiting integral is trivial to calculate, since Sp^p{x) — x, 
and so ^ 

lim riN{s) = ^5 

Thus, this sum shouldn't be interesting, except as a calibration of sorts. 

But that is false. A review of the figure [6] shows that this simple stair-step 
sum retains traces of the locations of the non-trivial Riemann zeros! 

The presence of this shadow is not hard to explain. Define 

^ 1 

n— 1 

Then it is clear that, for large A^, Zjv(s)will approximate the Riemann zeta 
function: 

lim Zt,{s)^C{a) 

A/— >oo 

at least for the case of 5Rs > 1, where the sum is convergent. With some simple 
manipulations, one may rewrite the numeric sum as 



/ 2n + 1 \ 1 



^ \ IN ) (2Ny 

n=0 ^ / \ ) 



[Z2N-l{~s)~rZN-l{-s)] 
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Figure 2: Riemann Sums 




The above strip charts show \Cn{s)\ for N=15, 151, 1051, 11051 and 151051 
respectively. The three vertical black lines are located at 5is=0, 1/2 and 1. 
The vertical range of these strips are all identical, ranging from 3s=13 at the 
bottom, to 5s=34 at the top. The coloration is such that red represents values 
greater than two, green for values about equal to one, and blue/black for values 
near/at zero. In this particular strip, we expect the exact Riemann zeta to 
have five zeros, near 14.13, 21.02, 25.01, 30.42, and 32.93. Several remarkable 
properties are visible: First, there are far more than five zeros visible, some lying 
outside the critical strip. Next, as the value of N increases, the zeros migrate 
towards the 5is = 1/2 line; however, the number of zeros seems to multiply 
logarithmically as well. Most of these zeroes seem to have a tiny residue, which 
diminishes as N increases; presumably, they will completely disappear in the 
limit. What's left are five zeros which are not fading away; these become the 
Riemann zeroes. 
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Figure 3: Riemann sum phase plots 




The above three figures are parametric phase plots of Cn{s) for N=1051, 11051 
and 121051, treating s as a parameter, varying from 3s=13 to 34 along the 
3?s = 1/2 line. That is, the graphs show on the x-axis, and Q(n along 
the y-axis as s varies. The large number of zeros in the critical strip can be 
understood as being due to the looping seen in these graphs. As N increases, 
the loops shrink, passing through cusps; the last figure beginning to resemble 
the true phase portrait of the Riemann zeta. 
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Figure 4: Zeros of Si^2 




This figure shows a side-by-side comparison of ICAfl^)! and |C-/V;i,2(s)|, with the 
later as defined in equatiorlSl and N — 151051. The sizes, heights and color 
scheme are as in the previous figures. Note the near-correspondence of the 
dominant zeros; yet, the zeros for |Cjv;1,^s)| do appear to be at shifted locations. 



Figure 5: Zeros of 81^3 and 5*2,3 
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This figure sliows a side-by-side comparison of |Cjv;1.3(s)| and |CAf;2,3(s)|, with 
the later as defined in equatior(3l and N = 151051. The sizes, heights and color 
scheme are as in the previous figures. The dominant zeros no longer resemble 
those of |Civ(s)|, and it is no longer cleag whether they will end up in the critical 
strip. However, it is also the case that the structure of Si^3{x) and 5*2, 3(0;) is 
much finer and detailed, and so this value of N does not yet provide an adequate 
visual estimate of the N 00 form of these figures. 



Figure 6: A Shadow of Zeta 




This figure compares |Civ(s)| and \riN{s)\ for N = 1105f. Note the "horns" or 
"cusps" in \r]N{s)\. While these are fahing well outside the critical strip, it is 
particularly remarkable that they line up vertically with the locations of the 
non-trivial Riemann zeros. ^ n 



The fact that the seemingly "trivial" rjN^s) sum contains this approximation to 
the zeta would seem to explain the numeric features shown in figure [6l Fur- 
thermore, the presence of the — s in the arguments appears to explain why the 
horns are located not far from the 5Rs = —1/2 line. 

On the other hand, there is another important lesson to be drawn here: this 
"trivial" sum is insufficient to explain the features seen in the CiV;p,g(s) sums. 
That is, although the numerical evaluation itself is prone to introducing the 
Riemann zeta as an artifact, the strength of this artifact does not immediately 
appear to be dominant enough to account for the entire behavior of the CN;p,q{s) 
sums. 



5 Sums and Integrals 

This section explores analytic approaches to evaluating the integral [51 Because 
Sp^q (x) is a piecewise combination of Mobius functions 

ax + b 



cx + d 

it is enough to evaluate the indefinite integral 

ax + b 



x" dx 



This integral does not have any simple, closed-form solution, but can be given 
in terms of an infinite sum. This may be obtained via the Newton series 



where 



fe=0 

Tis + 1) 



r(fc + i)r(s-fc + i) 

is the binomial coefficient. The series is absolutely convergent on the unit in- 
terval < X < 1 that we wish to explore. The indefinite integral may then be 
evaluated as a sum: 



F (s, a; x) 



dx 

X + a 

{y - aY 
y 



dy 



fc=0 



— ) ly'~'dy 



i-aT 



\n{x + a) 



k=l 
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This allows the Mobius integral to be written as 

ax + b 



I 



CX + d 



'-'^dx= -F {s,-;x] +-F {s-l,-;x 



Using that fact that, in general, for continued fractions, the a, b, c and d will be 
integers, with ad — be = ±1, one may re-write the above as 



fe=i 

(4) 



For 51^2 (x), this can be made more explicit. So, we have that 

'S'1,2 (a;) = / 1 

a2 + 



oi+r 

when X is written as 

1 

X = T— 

with ai = [l/xj and 02 = [x/ (1 — aix)\ integers. Solving for r and replacing, 
one gets 

ax + b 

Si, 2 {x) = — 

CX + d 

with 

a =1 + ai (a2 — 0.1) 
b =ai — a2 

c = {a2 - fli) (1 + 0102) 
d =1 + a2 (fli - 0,2) 

Note that in this case, ad — be = 1 for all possible values of ai and 02. This 
form is valid over the interval 



(x, x) = 



a2 l+a2 



_l + aia2 l + ai+aia2 
and so one may combine these elements to write 

i^"" ax + b 



/s,, ^..-.^ 

ai — 1 a2 — 1 — 



(5) 



At this point, the ability to perform any further simple manipulations stops. 
The summations are only conditionally convergent; one cannot interchange the 
sum over k with the sums over 04, 02. If one were to do so, one would quickly 
discover that, for fixed k, the sums over ai, 02 are logarithmically divergent. 
While this is not immediately obvious just be gazing at the sums, it is easily 
confirmed by numerically evaluating them. 
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One more avenue suggests itself, but then fails: one might consider applying 
the binomial expansion to the powers (1 + ^) that appear in the sum for 
F{s,a;x), and then exchanging the order of summations there. But this is 
not possible; the sums fail outright. Thus, it appears that there are no further 
(obvious) analytic techniques that one can apply to reduce the expression [SJ 
One is then left to contemplate the numeric evaluation of [5l Naively, one might 
expect that doing so would offer considerable advantages over the numeric sum 
[3l but this is not at all the case; convergence is even slower, and fails more 
profoundly. 

Thus, the remainder of this text uses only the sum[3]for numeric exploration. 
The position of the first non-trivial zero of Ci,2 (s) is examined in the next 
section. 

6 Locating the first zero 

This section explores the numeric evidence for the location of the first zero of 
Ci.2 (s)- One may quickly discover that it must be near s — ^ + i 14.92, but 
improving upon this value is remarkably difficult, given only the tools [3] and [S] 
The summations |3] are straightforward to evaluate numerically. Estimates 
for the location of the zero may be obtained using more-or-less standard zero- 
finding techniques, such as Powell's algorithmfS] (some care must be taken, 
as the neighborhood of the zero appears to be a bit wobbly). The resulting 
estimates for the zero depend very strongly on iV, and are very noisy, as the 
figure [7] shows. The noise appears to be classic fractal 1// noise, as can be seen 
from the graph of the power spectrum. 

7 Conclusions 

Numerical evidence suggests that the permutation-inspired zeta 

Cp,q{s) = lim CN:p.q{s) 

obeys a Riemann hypothesis, in that it's zeros lie in the critical strip, and partic- 
ularly on the IRs = 1/2 line. High-precision numerical computations are difficult, 
and offer only ambiguous support for the conjecture. Curiously, the numerical 
evaluation itself "accidentally" introduces some zeta-like artifacts, these artifacts 
do not seem to be sufficiently strong to explain the overall structure. Alternate 
approaches to the problem are desired, but the author is not aware of any par- 
ticular theoretical framework that would be applicable, and could be employed 
to provide deeper insight into the phenomenon. 
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Figure 7: Numerical zero-finding 
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The above graphs show the location of the first non-trivial zero of Cn:1,2 (s), as a 
function of N, in the range of 90000 < N < 95000. As may be seen, the location 
of the zero depends very strongly on N, and is an extremely noisy function of 
N. There is the vaguest hint of an oscillatory behavior in N, possibly with a 
period of AA^ = 300, but this is very strongly obscured by the noise. The power 
spectrum of this same data-set is examined in the next figure. 
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Figure 8: Power spectrum 




These two charts show the power spectrum of the zero locations graphed in 
figure [71 The power spectrum is defined as usual, it is the square of the Fourier 
transform of the data. The straight line shows the 1// spectrum; the noise in 
the location of the zero clearly obeys the classic fractal 1// power density. The 
regular pattern of arches appearing at low frequencies is a numerical artifact 
resulting from the limited size of the data sample. The left chart shows the 
power along the real axis (resulting in the uncertainty of the 5Rs = i location) , 
while the right chart shows the power in the imaginary direction. The higher 
power in the right chart indicates that the location 5s w 14.92 is more certain 
than the location 5fts « i. Curiously, the "more important" location is harder 
to pin down! 
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Figure 9: Increased precision 



Zero of ?etH_N;1 ,2 Zero of ?etH_N;1 ,2 




The above graphs show the location of the first non-trivial zero of Ca';1,2 (s), 
as a function of N, on a logarithmic scale. Only a sparse set of values of 
N are sampled. The figure on the left, which would confirm or disprove the 
conjecture, is ambiguous. Although the data, as illustrated, shows a precipitous 
drop past the critical line 5Rs = i, it is not yet obvious whether there will 
be a bounce/oscillation back up towards 5Rs = i, or whether the asymptotic 
behavior aims at 5Rs ~ 0.45. Unfortunately, computational costs are large; each 
additional point in these graphs requires in excess of a week of compute time 
on current-era computers. 



[2] A. Ya. Khinchin. Continued Fractions. Dover Publications, (reproduction of 
1964 english translation of the original 1935 russian edition) edition, 1997. 

[3] Wilham H. Press, Saul A. Teukolsky, WiUiam T. Vetterling, and Brian P. 
Flannery. Numerical Recipes in C, The Art of Scientific Computing. Cam- 
bidge University Press, 2nd edition edition, 1988. 

[4] Linas Vepstas. On the minkowski measure. ArXiv, arXiv:0810.1265, 2008. 



16 



